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Abstract 

We construct a phenomenological theory of gravitation based on a second order gauge formulation 
for the Lorentz group. The model presents a long-range modification for the gravitational field 
leading to a a cosmological model provided with an accelerated expansion at recent times. We 
fit the model parameters - such as the coupling constant of Lagrangian's sector scaling with the 
square of V-R - using super-novae type la and X-ray galaxy clusters data. There is good agreement 
between our estimative for the age of the universe and the one predicted by the standard model. 
The transition from the decelerated expansion regime to the accelerated one occurs recently enough 
to indicate a solution of the cosmic coincidence problem. 



I. INTRODUCTION 



To explain the origin and evolution of the present accelerated expansion of the universe is 
one of the most challenging problems of Physics nowadays. One way to obtain a mechanism 
of acceleration is changing one of the cornerstones of modern physics, the theory of General 
Relativity. 

Modifications in the scheme of General Relativity are being proposed since its invention, 
in the beginning of the 20th century, and they are motivated by several reasons, from the 
quest to agreement with the theory for the inner structure of quantized matter to eventual 
need of extra-dimensions and the desire to achieve the unification of interactions. The first 
modification of General Relativity was proposed by Einstein through the introduction of 
the cosmo logical constant, which is one of the several alternatives to describe the present 
acceleration of the universe. Other proposals in the context of Cosmology involve the in- 
troduction of one (or more) spatial extra-dimension in the braneworld scenario p], or the 
presence of a self- interacting scalar field, the quintessence models [2]. 

Here, we shall explore some of the cosmological consequences of a phenomeno logical the- 
ory of gravitation based on a Lagrangian analyzed in [3]; we will show that this theory 
permits a recent accelerated phase for the universe without the introduction of A, exotic 
matter, extra-dimensions or scalar fields. It was constructed on the basis of a gauge formu- 
lation for the gravitational field, through the second order gauge theory [Ij. This theory, 
by its turn, was shown to be efficient in the description of effective limits of other gauge 
theories, as Podolsky electrodynamics and the SU{N) non-abelian model. 

Our model uses differential equations of order higher than two in the metric tensor, but 
it is not equivalent to the usual approaches quadratic in the curvature tensor [5]. It is also 
not equivalent to the f{R) theories discussed, for instance, in |6] and [7]. From the gauge 
theoretical point of view, quadratic Lagrangians in the Riemann tensor are actually of first 
order, and a classification of all possible quadratic Lagrangians of first and second order 
in the gauge gravitational field are given in [3]. Among these possibilities, we choose a 
particularly simple one, which is inspired by the Podolsky's abelian case and by the effective 
Alekseev-Arbuzov-Baikov's non-abelian model [8] . Actually, we choose to add a term scaling 
with the square of the covariant derivative of the scalar curvature V-R to the familiar R of 
the Einstein-Hilbert Lagrangian. Because of this choice, our study is not the most general 
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one: we do not take a linear combination of all the gauge invariant Lagrangians presented 
in [3]. Our main intention is to analyze the consequences of this option in the cosmological 
context, and in particular if higher order derivative terms could produce interesting effects 
- such as the present-day cosmic acceleration of the scale factor. 

We will adopt the following action for the description of the gravitational field: 

S = I d'x^g (^^ + ^Cp - , (1) 

Cp = Iv.RV^R . (2) 
o 

Notice our option of not including the cosmological constant term in the action. At this 
stage this is only a simplicity hypothesis, but - we repeat - our final goal is to reproduce the 
present-day cosmic acceleration effect produced by A without imposing its existence a priori. 
The final results will show if our choice is paid of. It is worth to emphasize that the choice 
(|2| constitutes a phenomenological model valid within a limited interval of energy (set by 
the values of the coupling constant it does not hold during all the cosmological history 
(as we shall see) but only for a certain period. The same phenomenological Lagrangian was 
applied to describe inflation in Ref. [9], but we emphasize the high order of magnitude of 
the energies involved there, which would correspond to the early stages of evolution of the 
Universe. On the other hand, the model presented here engenders sensible dynamical effects 
at recent cosmic times, where the energy scale is very low. 

The paper is organized as follows. In section ITT] the field equations are written for a 



Friedmann-Lemaitre- Robertson- Walker metric. Section III is devoted to obtain a perturba- 
tive solution of the field equations about the usual dust-matter model of the Einstein-Hilbert 
theory (as described by the Friedmann equations). The perturbative solution is constructed 
in such a way that the universe is dominated by a decelerated regime until the time t* from 



which the additional term begins to be relevant. In section IV the parameters of the 



model are related to the observational data, through a set of coupled nonlinear equations. 
Such equations are solved by numerical methods in section [V], and the results are discussed 
in section |VT1 
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II. FRIEDMANN EQUATIONS 



The invariance of the action ([T]) with respect to 5gxu yields the field equations: 



Hxu = VaV. [OR] + ^VxR^uR - RxuOR - gxuO [OR] - ^gxu^'RVpR • 



(3) 



where = V^V^ and is the covariant derivative. 

Applying the field equations (jsj) to a homogeneous and isotropic space, described by the 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, 



1 



-dr' + r'^dO? 



1 — /tr^ 

(k = — 1,0,+1) one finds, after some direct but long calculations 



-3HR + 3HR - GH^ii + 9H^R + \r^ 

4 



00 ; 



'1 — Kr^ 



H + 3H^ + -R + 2— + p("R+ 5H'r + 3H^R + 5HR+ 



-3HHR - 9H'^R + 3HR +-R^-2^{R + 3HR 

4 
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(4) 



where H (t) = d/a is the Hubble function, R{t) = g^^R^u is the scalar curvature and we 
are using units such that x = SttG. In our notation, dot means derivation with respect to 
the cosmic time t. These are the higher order Friedmann equations in terms of the Hubble 
function H (t) and the scalar curvature R (t). 

Following the standard procedure we use the energy-momentum tensor of a perfect fiuid 
in a commoving coordinate system, 

Tf,^ = {p + p) ^Jf^v - P9^.u ■ 

In order to simplify the treatment, we will be concerned only with the case of a fiat 
spatial section, k = 0. So, using the relationship between the scalar curvature and the 
Hubble function. 



R 
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we get the following modified Friedmann equations: 

3//2 + p (18HH + lOSH^H - 18HH + 9H^ + 90H^H- 

+216HHH - 72H^ + 288 (m)^ - 216H^H^ = xp , 

2H + 3H^ + (3 (gH^^^ + 5AHH + ISSH'^H + 12QHH+ 
+81H^ + im^H + A98HHH + 120H^ - 216H^H 



-XP ■ 

Combining the equations, one finds: 

2H + (3 (qh'-^^ + 3QHH + 30H^H + lUHH + 72H^+ (5) 

-12H'^H + 282HHH + 192if^ - 288 (m)^^ =-x{p + p) ■ 

As we want to describe the evolution of the universe, this equation must be complemented 
with the covariant conservation of energy-momentum, 

p + 3Hip + p)=0 , 

and an equation of state / relating the energy density p, the pressure p and the Hubble 
function H, 

fip,p,H) = 0. 

The dependence on H is included to account for the general case when one admits interaction 
among the constituents of the cosmic fluid [H]. In this case, there is a possible constraint 
relating p, p and the scale factor, or equivalently H. On the other hand, the usual equations 
of state of physical cosmology associate only pressure p to the energy density p, or pressure 
to the numerical density n. For example, the equation of state for the dust matter is 
p = nkT -C p (A; is the Boltzmann constant and T the temperature), and p = p/3 is the one 
used for ultra-relativistic particles. 

III. SOLUTIONS OF THE HIGHER ORDER FRIEDMANN EQUATIONS 
A. Dust Matter 



Our main interest here is to apply the model to the present state of the universe. There- 
fore, we take as source a perfect fluid composed by dust matter p = (ordinary or dark). 
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In this case, the continuity equation gives: 

^^')=^°(^) • 

In order to use such result directly, we would have to rewrite equation ([s]) in terms of 
the scale factor, obtaining a nonlinear and much more complicate equation, which we shall 
avoid. Instead, we consider simultaneously the following pair of coupled equations: 

H + 13 (3H^^^ + 18HH + 15H'^H + 72HH + 36H'^+ 

-36H^H + UIHHH + 96H^ - 144 (m)^^ " "f*^ ' '^^^ 

p+3Hp=0 . 

These equations can be analyzed by several methods, such as linearization for dynamical 
systems, spectral analysis or perturbation theory. Here, we will consider only this last 
procedure, leaving the other options for future investigations. 



B. Perturbation Theory 

The model is constructed by assuming a standard Friedmann expansion prior to some 
time t* from which the second order effects start to become significant. The strategy is 
to consider a perturbation series in the coupling parameter j3, in order to guarantee the 
accordance of our model with the usual cosmological model (in some region of the space of 
parameters). Take, for instance, an expansion up to second order terms; it reads: 

H (t) = Hf + pHi + (3^H2 , 

p (t) =PF + PPI + P^P2 , (7) 

where the label F stands for the standard Friedmann solution of the Einstein equations. 

Substituting expansions ([t]) in the pair ^ and matching the terms order by order, we 
get: 

0(,3»)^|^-+i''- = ° . (8) 

Pp + 3HpPf = 



and. 



0(^.)_.^"' + i'>' = ^'W , (9) 

Pi + SHfPi + 3HiPf = 



where 



and also, 



with, 



Si {t) = -{ 3HP + ISHfHf + I^HI'Hf + 72HfHf + SGH^f 



-36HiHF + UIHfHfHf + 96M - IUHIhI 



P2 + 3HfP2 + 3H2PF = -^Hipi 



(10) 



(11) 



^2 (t) = - ^Hf^ + 18 {HiHf + HfHi) + 30HfHiHf+ 



72 HfHi + HiHf + HfHi - IOSHiH^Hf - 36H^fH 



(12) 



+ 

141 ^HiHfHf + H fH fH\ + HfHiHf j + 288 ^HiHp — H^H fH\ — HfHiHf 

This way, we obtained a pair of coupled linear equations in each order. Their previous 
orders give the source term and the coefficients. 

1. Zeroth order solution: the standard cosmological model 



The solution for the system of zeroth order in the coupling parameter /?, Eq. ([8 

i^F + |PF = , (13) 

PF + 3HfPf = , (14) 



can be obtained by direct integration. Substituting (13) in (14), we have: 

3 2 J 
Pf = -Hp + ki 
X 



where ki is an integration constant. Substituting this back into (13) and choosing the initial 
condition /ci = 0, 

Hf 2 ^ 

Setting = 0, we obtain the standard Friedmann solution for H {t): 

2 1 



3t 



(15) 
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2. First order solution 



In the first order approximation, we have the coupled set ([9]), 

Hi + |pi = S, (t) , 

pi + SHfPi + SHiPf = , 

These equations can be solved by the Increasing Order Method. Differentiating the first of 
these equations and using the second one, we obtain 

Hi + 3HfHi - yPfHi = Si (t) 
Si (t) = Si (t) + SHpSi (t) 

The general solution of such equation can be obtained in the form of a power law: 

rr / N ,2 4912 . 
Hi (t) =at + bt~^ H 

, , 2 / 15 977 , o\ 

The integration constants a and b should be chosen in accordance to the physical situation 
to be described. Since the zeroth order terms appear as source terms in the first order 
approximation, one can choose the initial conditions Hi (t*) = pi (t*) = 0. This determine 
the integration constants leaving the theory with only three free parameters, namely the 
coupling constant /5, the age of the universe to (see below) and the instant of perturbation 
t*. 

Therefore, in the first order approximation we find the following solution to equations 
Q: 

H(t)=--+^( ^-^ t-' + — f-Yt-' - — - (t*)A 

^ ' 3t (t*)^ \^ 729 V ^ y 243 V V 729 t* ^ ' J 

8,GpW = ii + 2A('— f-V-^- — f-V*-^ + — (O-^) (16) 

'^^^ 3t^ (t*) V 729 V^/ 243 V^/ 729 ^ ^ ^ ^ ^ 

IV. OBSERVATIONAL PARAMETERS 

Now, let us focus on the problem of linking our theoretical model with the observational 
data available. Once the main goal of the present model is to reproduce the present accel- 
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eration without using the cosmological constant or scalar fields, we must choose the data 
which are model-independent. 

The redshift z and the luminosity distance di^ are dependent on the null geodesic equation 
only, and this is not changed by the second order field equations ([s]). Therefore, they con- 
stitute the ideal data set to be compared with the predictions of our model. The luminosity 
distance can be directly related to the redshift [T2] . 

1 ( 1 



~ + ^ (1 - go) 

go is the deceleration parameter. 

The supernovae projects usually measure the curve of rf/, (z) determining the parameters 
i^o and go with good accuracy. However, in our model we have three parameters to be found: 

1. The age of the universe to; 

2. The instant t* from which the perturbation coming from the modified gravitational 
equation becomes important; 

3. The coupling constant /? for the higher derivative terms in the action. 

We need three independent measurements to find these parameters. Two of them are 
naturally given by the supernovae experiments ~ through the values of i^o and go- The third 
is suggested by the very conservation equation for the energy momentum tensor. Since the 
presented model considers only dust matter, we need the value of the total amount of matter 
in the universe, or equivalently its ratio with respect to the usual critical density, fimo- This 
information can be provided by X-ray experiments measuring the gas fraction in galaxy 
clusters [T3|. The advantage of using the results of this particular measurement is the fact 
that it is completely independent of the model for the evolution of the universe. This is not 
the case for the more common WMAP data: they assume the validity of h^CDM model. 

In the following section, we shall carefully discuss how to use Hq, and to obtain 



to, t* and p. But, before that, we will add to system (16) the constraint 

Hito) = -H^iqo + l) (17) 
following from the definition of both the Hubble and the deceleration functions in terms of 



the scale factor: H = a/ a, q = —da/a^. Gathering (16) and (17), we get the new system to 



be solved: 



2 1 ( 35408to.,..-i 11065 /to V\-i^ 4912 ^to 



3 to {t* f\ 729 t* ^ ^ 729 V^V 243 

• 2 1 P ( 35 408 2 22130 /to 

S^r.rM ^^^9^^ f 15977 /toy\_2, 35408 2 ,22130 /to 
87rGp(to) = -;2+2— ^ -— — - + 




3t^ (t*)M 243 V^*/ 729 ' ' 729 

The first member of each equation of the system above is given in terms of observational 
constants, while the right hand side of each equality bears the parameters of the perturbed 
model. These will be calculated by solving numerically the above transcendental equations. 
We deal with this task now. 



V. NUMERICAL CALCULATIONS 



In order to solve numerically the system of coupled transcendental equations, let us 
perform some simple manipulations. First, we define new non-dimensional variables, 

t* B 

u= - ,b=^; (18) 
to ' t*4 ' ^ ^ 

in terms of which the system is rewritten as 

2 1/. , / 17704 1 11065 2456 4 



^»=3fcil + H-Wi^ + ^"+^"JJ^ '''' 
X 2 1 / /17704 1 11065 12280 A\ 

. 4 1 / / 15 977 4 17704 1 11065 



where 



^mO = PO ■ 



Taking the ratio of the two last equations of the system above, one gets: 



2jqo + ll^Sn^^ (22) 
s (u) 



, , , /17704 1 11065 36840 A 40537 , , . 

. (u) ^ (3O.0 - 2 (,o + 1)) ^ + + —u^J + ^ (^0 + 1) . 
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We must have P < in order to assure the stabihty of the theory jSj |9]. This estabhshes a 
constraint on the sign of b. Substituting b {u) in the equation for Hq: 

1 3 Ho 



toiu) 2il + biu){--^± + -^u + Wu^)) 



(23) 



Combining (19) and (20) we get 



r2.o.„ , 1^ on , _ 4 & (^) /12280 15977\ , 

3tl{u) V 81 162 J ' 



H^i2 (go + l)-3fi^o) 



which is a nonhnear equation for the parameter u. This equation can be solved using the 



Newton- Raphson method. The result is then used in (22) to obtain the numerical value 



of b and both b and u are inserted in (23) to get the age of the universe to- The time 
of perturbation t* can then be obtained as a simple ratio. Proceeding this way, we are 
able to find all the parameters of the model given the measurements of Hq, go and po (or 
flmo)- Notwithstanding, such measurements have uncertainties, which are propagated to the 
parameters of our model with a confidence level of 68% (one a). This involves a problem 
of nonlinear programming [T7| using the algorithm described above implemented in the 
software Mathematica® 6.0. 

With this assumptions and using the experimental values [131 E] 

Qo = — 0.81^0^4 , 
O — n 'ici+o oe 

■^'mO — '-'•'-''-'-0.06 ' 

we find: 

6 = -0.0037l°:r, 
t* = 9.A+\i Gyr , 

which can be used to calculate the parameter j3 and the age of the universe, 

P = -29tZ {Gyrf , (24) 
to = Gyr . (25) 

In spite of its high value, [3 does not break the meaning of the modified action as proposed 
in ([T]). This is so because the own curvature tensor (and its derivatives) is /3-dependent, such 
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that the additional term Lp = ^^/—g {f3)'V^R V^R {(3) scales as and therefore is 

small when compared to the usual Einstein-Hilbert term. The fact that the uncertainty 
in (3 is very large - it is of the same order of /5 - is not critical for our model precisely 
because <^ 1. Of course, the contribution of the additional term goes up for values 

of t progressively greater than t*. But we do not hope that the perturbative approach 
holds for arbitrary large values of t. The convergence of the perturbative expansion can 
be qualitatively studied plotting the ratio of the perturbation term by the usual FLRW 
solution, as given in Fig. [T] 




t(Gyr) 



FIG. 1: Ratios of the first term in the perturbative solutions by the ordinary dust-matter solution as 
functions of the time. The fuh hne shows the evolution of the the ratio of the first-order perturbed 
Hubble function divided by the non-perturbed one; the dashed line presents the behavior of the 
perturbation on the energy density divided by its non-perturbed values. 

As long as the curves lie within the interval [—1, 1] the perturbative scheme can be 
assumed as valid. Notice that the Hubble function fits the perturbative scheme for a longer 
time compared to the behavior of the energy density. 

The value of to provided by the WMAP3 data [16] - which assumes the ACDM model 
- is calculated as 13.7]'^o2 Gyr (for a fiat Universe). The result given by our model is 
in agreement with the one predicted by the ACDM model at the 2(T-level. This apparent 
tension between both results does not impair the model proposed here since results obtained 
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from anisotropies are strongly dependent of the adjusted model, in this case the ACDM one. 
Besides, globular clusters data gives 11.2 Gyr, at 95% of confidence level, as inferior limit 
for the age of the universe [T5] . 



A. The "Coincidence Problem" 



Using the solution (16), one easily finds the ratio of scale factor at two arbitrary times ti 
and t/ as 



ln^ = ln^^^ 



13 { 17704 ftfY ^ 11065 ftf 



1228 /t,^ 



Uj (t*)^ \ 729 \t* J 729 
p ( 17704 fuY^ 11065 fUy^ ^ 1228 f U 



243 \t* 

4' 



+ 



{t*f \ 729 \t* ) ' 729 \t* J ' 243 \t* 
which exhibits a very smooth transition from the Friedmann standard regime to the accel- 
erated one - see Figure bellow. 




t (Gyr) 



FIG. 2: Plot of the complete scale factor (matter solution plus higher order terms) as a function 
of time. Notice that the decelerated regime is roughly in the first half of the graphic; the curve 
shows clearly the onset of acceleration at recent times. 



The same quantity can be used to estimate the red-shift at the transition, 

ao 



1 + z* 



a{t*) 
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as 

(to\ \ ^ ( P ( 17704 /toy 11065 /tpV^ 1228 /tpy' 32A53 \\ 

^ \t* J ^^^^ 1^ (t*)^ y 729 V^* J 729 V^* J 243 \t* J 729 J J 

With the values of /3, t* and to, one calculates: 

^ — U.iJ_o.07 ■ 

Therefore, our second order Lagrangian-based model presents a possible explanation for the 
"cosmic coincidence" problem, since the observed values of the present rate of expansion 
{Hq), acceleration (go) and total mass (i^mo ) imply a transition from a decelerated phase 
to the present accelerated regime at a relatively low redshift. 



VI. CONCLUSIONS 



We have constructed a model based on a phenomenological theory of gravitation obtained 
from the inclusion of a Podolsky-like term scaling with the square of the covariant deriva- 
tive of the Ricci scalar. This model implies long-range modifications in gravitation, which 
leads to an accelerated regime for the present-day universe, even in the absence of a dark 
energy component or cosmological constant. According to our perturbative evaluation, this 
accelerated expansion started recently as indicated by the values of t* or, equivalently, z*. 

The values predicted by the model for the age of the universe and the redshift of transition 
are in agreement with the supernovae data ^Mj or the analysis of the cosmic microwave 
background based on the ACDM model at the 2a-leve\. We did not use the CMB data 
to predict values for the age of the universe, our intention being to avoid model-dependence. 
But these calculations could be done in the future. Our perspectives include the study of a 
perturbative solution for a closed {k = 1) universe, keeping in mind the exact solution found 
in the ordinary FLRW case [19|. 

The fact that the gravitational field is weaker at long distances with a characteristic scale 
given by the coupling constant (3 suggests the existence of massive modes in the weak field 
approximation, but in a way that does not break the coordinate invariance, analogously to 
what happens in the Podolsky electrodynamics The eventual existence of such massive 
modes are under investigation, and the results should be compared to other approaches in 
the same direction l20l. 
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